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Abstract 



We overview high-energy neutrinos from galactic sources, transparent to their gamma-ray emission. We focus on young supernova 
remnants and in particular on RX J1713. 7-3946, discussing expectations and upper bounds. We also consider the possibility to 



detect neutrinos from other strong galactic gamma-ray sources as Vela Junior, the Cygnus Region and the recently discovered 
Fermi Bubbles. We quantify the impact of the recent hint for a large value of #13 on high-energy neutrino oscillations. 
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C3 1. Context, motivations and assumptions 



The successes of low energy neutrino astronomy and the 
discovery of neutrino oscillations added momentum to the 
search for high energy neutrinos from cosmic sources, initi- 
ated long ago with the theoretical proposals of Zheleznykh, 
Markov and Greisen [1]; see Fig. 1. The first km 3 -class 
detector, IceCUBE and the smaller ANTARES (follow- 
ing MACRO, BAIKAL, AMANDA) have not yet revealed 
the first signal, yet the excitement remains high. However, 
nowadays it becomes clear that the search for high energy 
neutrino sources is difficult. In this paper, we attempt to 
examine certain generic and specific expectations for such 
a search, in the hope to contribute to the scientific plan- 
ning of the future experimental activities. We are inter- 
ested in the possibility to proceeding further, by building a 
new telescope of km 3 -class in the Northern Hemisphere-as 
Km3NET in the Mediterranean Sea or GVD in Lake Baikal. 

The fact that the first explorations of the high-energy 
neutrino sky, that did not reveal any signal, contradicted 
the over-optimistic expectations may lead one to doubt that 
we can reach reliable predictions. But, even if some degree 
of diffidence toward theory is healthy, we are convinced that 
the process of formulating expectations is an important step 
toward understanding. We would like to clarify our view by 
proposing some issues: Do we really need to have high expec- 
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tations for high energy neutrinos? Perhaps not, although 
one should keep in mind that the surprises often happen 
in astronomy. As recent example in this regard are the dis- 
covery of the Fermi Bubbles and the variability of the Crab 
Nebula. Are expectations useful? Of course, yes; good pre- 
dictions are precious for experiments, but also reasonable 
expectations eventually contradicted are not useless. Do 
we have some relevant and succesful precedent? Yes, many: 
High-energy 7-ray sources have been predicted in the fifties 
(Morrison). For solar neutrinos, we had predictions with 
errorbars since the sixties (Bahcall) . For supernova neutri- 
nos, quantitative expectations were available even before 
SN1987A (Nadyozhin). 

We adopt a rather common astrophysical attitude: the 
observation of high-energy neutrinos is very demanding but 
perhaps possible and would amount to an unambiguous sig- 
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Fig. 1. Energy ranges of various types of neutrino telescopes, in 
log 10 (E /eV) ; note that they span more than 10 decades in energy. In 
this paper, we are concerned with the third item, that after IceCUBE 
we identify with the technology 'ice Cherenkov'. 
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Fig. 2. We can exploit the strict connection of secondary muon 
neutrinos and ■y-rays produced by collisions of cosmic rays above 10 
TeV, assuming that the ■y-rays are not absorbed or strongly modified. 

nal of cosmic ray collisions-hopefully, those in their source. 
We simplify our task by restricting our attention on the 
specific but important class of sources, those transparent 
to their gamma rays. We will show how to proceed towards 
precise expectations and eventually observations of neutri- 
nos, limiting the use of theory inputs and just by the help 
of 7-ray observations. This aspect is important enough to 
be emphasized, as we do with Fig. 2. In the next section, 
we quantify the connection between gamma and neutrinos. 

2. Expected neutrino intensity and signal 

If we measure the very high 7-ray emission from a cos- 
mic source, and if we attribute it to cosmic ray colliding 
with other hadrons, it is straightforward to derive the muon 
neutrino flux. In fact, when both the neutrinos and unmod- 
ified, hadronic 7-rays are linear functions of the cosmic ray 
intensity, they are linked by a linear relation [2] : 

M E)= 0.380 /,( I ^r)+ 0.013 Lr(^) (D 
+ / f (!) 



where the coefficients are determined by the amount of 
mesons produced and by their decay kinematics. The first 
and second contributions are due to direct decays into neu- 
trinos and the third one to [i decay, the kernel K M (x) being 

' x 2 (15.34- 28.93x) < x < r K 

< 0.0165 + 0.1193x + 3. 747x 2 - 3.981x 3 r K <x<r w (2) 
^ (1 - x) 2 (-0.6698 + 6.588x) rv < x < 1 

where r x = {m^/rrix) 2 with x — ir,K. Similar formulas 
describe antineutrino production; the effects of the oscil- 
lations are included [2], compare also with [3] and see the 
Appendix for an updated discussion. 

At this point, calculating the muon signal is just applying 
a standard procedure. First, one evaluates the probability 
of converting muon neutrinos into muons 

E 

p v^p.= J dE^-^R^/m n (3) 
[say, 10~ 35 cm 2 x /m„/3 ~ 10~ 6 ] 



dN,, 

E * [A.U.] 

dE 




Fig. 3. Distribution of leading to muons, assuming E~ 2 pri- 
mary spectrum (upper curve); then, including Earth absorption, for 
a source at 5 = —39° as seen from ANTARES (middle curve); then 
with a spectrum E~ 2 exp( —yj E/150 TcV ) (lower curve), i.e., with 
primaries cutoff ed at ~5 PeV. 

Unless stated otherwise, we will use E t h = 1 TeV in order 
to cope with atmospheric neutrino background; see [4] and 
Fig. 1. Then, one should calculate the effective neutrino 
detector area, proportional to the previous probability, to 
the physical area of the muon detector (that in general 
depends on the zenith angle 9), and to the probability that 
the Earth absorbs neutrinos: 

= A^O) x P V ^{E, 9) x e~° */™" (4) 
[say, 1 km 2 x 10~ 6 ~ 1 m 2 ] 

When E ~ 10 TeV, s - 2m n E - Q 2 > M w , thus this 
is the energy where the cross section begins to grow more 
slowly. The Earth absorption begins at E ~ few • 100 TeV, 
when a(E) — m„/(i? ffi p ffi ) ~ 5 • 10~ 34 cm 2 , see Fig. 3. 

3. Neutrinos and gamma rays 

If one assumes a power-law distribution with an exponen- 
tial cutoff for the primary cosmic rays, then the spectrum 
of 7-rays has the same power-law index but drops slower [5] 

Ly cx E~ a ■ eM-\/E 7 /E c ). (5) 

Typically, in 7-ray sources a varies between 1.8 — 2.2 and 
and E c varies from 1 TeV to 100 TeV. Assuming 7-ray 
transparency we can calculate the critical (minimal) inten- 
sity that corresponds 1 muon/km 2 year above 1 TeV; this 
is given in Fig. 4, where we consider a = 1.8, 2.2 and E c = 
1, 30, 1000 TeV. Interestingly, the integrated flux above 10 
TeV is almost the same in all these cases 

I y (> 10 TeV) = (1 - 2) x 10~ 13 /(cm 2 s) (6) 

Note that the present Chcrenkov telescopes are less sensi- 
tive above 10 TeV, due to the small collection area. Thus we 
ask how to detect these fluxes. In order to collect > 100 7's 
in a reasonable time, we need a km 2 area: 

Exposure = L 2 x T - 2 x km 2 x 10 h (7) 

e.g., a 10x10 Cherenkov telescopes array, or one dedicated 
EAS array. A large area 7 apparatus, such as CTA or even 
a custom instrument, could be then worthwhile even for 
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Fig. 4. ■y-ray integral intensities, corresponding to a signal of 1 
muon/km 2 yr and above a threshold E^, evaluated assuming that 
the sources are transparent to their gamma rays. See the text, and 
compare with [6], where we show the differential intensities instead. 

neutrino community, and should cost much less than a 
neutrino-telescope. 

Now that we have compared the energy-response of 
gamma- and neutrino-telescopes, we would like to empha- 
size the fact that 7 and views are complementary, as sug- 
gested by Fig. 5. The maximal complementarity would be 
obtained with two detectors in antipodal locations. Stated 
in more quantitative terms, a steady source at declination 
8 is seen from a detector at latitude <f> for a fraction of time: 

/ 7 = Re[cos _1 (— tan(Stan0)]/7r; (8) 

the fraction of time for neutrinos is just f v = 1 — / 7 . 
A concrete way to illustrate this feature is to note that a 
hypothetical (resp., 7) emission from Galactic Center, 
6 m —30°, is visible from North (resp., South) Pole. 

This argument can be generalized. Suppose that the 
galactic sources of TeV neutrinos trace the mass distribu- 
tion, and calculate the relevance of a detector as a function 
of its latitude. The matter of the Galaxy is mostly located 
in the region 6 < 0, i.e., below the celestial equator. Thus, 
a telescope at the latitude of NEMO has a priori 2.9 better 
chances (or 1.4, weighting the mass distribution with 1/r 2 ) 
to see galactic neutrino sources than IceCUBE as shown 
in Fig. 6 and argued in [6]. These remarks show why we 
think that a large, high-energy neutrino detectors in the 
Northern Hemisphere, not only would complement the Ice- 
CUBE field of view, but would also be important to probe 
the existence of hypothetical galactic neutrino sources. 

However these considerations, as general as they are, 
leave something to be desired. For instance, the use the av- 
erage galactic matter distribution to infer the source neu- 




Fig. 5. The main particularity of neutrino-telescopes is that they look 
downward; in fact, due to atmospheric fi background, from cosmic 
sources are preferentially detected from below. 
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Fig. 6. Continuous line: detector relevance, calculated assuming the 
galactic mass distribution. Dashed: same, weighted with 1/r 2 to ac- 
count for the fact that nearby sources are easier to detect. From [6]. 

trino distribution is doubtful, also because HESS sky-maps 
show only few intense 7-ray sources. In this regard it is 
more appropriate to study individual sources. 

In order to proceed, it is useful to formulate explicitly 
some pending question and doubt: (i) Is 7-transparency a 
reliable hypothesis? (ii) Are we sure of the 'point-source' 
hypothesis? Similar as asking whether <C 1° pointing is re- 
ally important for very- high-energy 7 and/or v telescopes? 
(iii) 7 above 10 TeV can help v astronomy, and this would 
be a natural direction of progress anyway, but do we have 
any guaranteed aim for such a search? (iv) Eventually, the 
true question is: How to separate leptonic from hadronic 
gamma's? Neutrino identify CR collisions, but is this the 
only way to proceed? 

All these questions suggest to define more precisely which 
(type of) source we want to investigate. 

4. Supernova remnants+molecular clouds 

The hypothesis that supernova remnants (SNR) are main 
sites of the acceleration of the galactic cosmic rays (CR) 
can and should to be tested. 2 As indicated in ref. [7], the 
associations between SNR and molecular clouds (MC) pro- 
vide the optimal conditions to produce intense neutrino- 
and 7-ray fluxes, that thereby can be used as a tool to test 
the SNR paradigm. 

The best hope to detect 'hadronic' neutrinos is connected 
then to young SNR, i.e., of age less than 2000 yr. In fact, 
these are expected to contain the highest energy cosmic 
rays, and can thus produce energetic neutrinos that stand 
over the background of atmospheric neutrinos. In other 
words, we are interested in special types of SNR's as poten- 

2 Let us recall the main steps that consolidated this hypothesis: 
1) Baade & Zwicky's propose that supernova explosions and CR 
are connected; 2) Fermi argues that the kinetic energy of gas can 
convert into CR acceleration via magneto- hydrodynamical processes; 
3) Ginzburg & Syrovatskii note that the kinetic energy injected by 
SNR is 1-2 orders of magnitudes greater than the galactic energy 
losses of CR; 4) several authors note that the shock wave present in 
astrophysical environments as SNR offers the conditions for efficient 
CR acceleration; 5) Bell & Lucek et al. remark that the CR lead to 
an amplification of the relevant magnetic fields in the region of the 
shock wave, as observed; 6) the present efforts are directed to develop 
a full non-linear theory of CR and SNR possibly including MC. 
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Expected 
signal 


1 sigma sta- 
tistical error 


Atmospheric 
background 


50 GeV 


5.7 


6% 


21 


200 GeV 


4.7 


7% 


7 


1 TeV 


2.4 


10% 


1 


5 TeV 


0.6 


30% 


0.1 


20 TeV 


0.1 


100% 


0.0 



Table 1 

Dependence on the threshold of the number of signal rauons per 
km? per year from RX J1713. 7-3946 and assuming the hadronic 
hypothesis. Also quoted the estimated error from HESS statistics and 
the estimated background. From [2]. 

tial neutrino sources: the closest and the youngest, among 
those associated to molecular clouds. It should however be 
observed that core collapse supernovae, being connected to 
short living stars, are pretty naturally associated with sites 
of intense stars formation-i.e., molecular clouds. 

Before analyzing a specific and important case of such 
a young SNR, let us discuss in some generality the pros 
& cons of the SNR+molecular clouds paradigm, outlined 
above. Some support comes from the detected GeV 7's from 
relatively old SNR, as W28 and W44. Moreover, gamma 
transparency usually holds. One can expect, a priori, that 
the closest SNR should be at about 1 kpc; in fact, we have 1 
new SN each many tens or years and the size of the Galaxy is 
several tens of kpc. Also, it is true that the CR acceleration 
above 100 TeV is still an open theoretical problem but, on 
the other hand, it is a problem that can be approached 
observationally measuring gammas above 10 TeV. 

In any case, we should expect that in order to describe 
any individual SNR we need theoretical modeling anyway 
(and the presence of MC will not make the matter simpler) 
plus as many multiwavelength observations as possible. 

5. The SNR RX J1713. 7-3946: a case study 

In this section, we will discuss the object RX J1713.7- 
3946, in particular, in connection with the possibility that 
this is an intense source of neutrinos above TeV. The na- 
ture of the object, its age and its distance are still debated 
issues; however, the identification of the location of this 
source of non-thermal X ray emission [9] with the position 
of a historical supernova [10] supports the idea that this 
is a SNR about 1600 years old. The distance is estimated 
around 1 kpc [11]; thus, it is pretty close. Moreover, the 
observations indicate an existence of dense condensations 
inside the shell of RX J1713.7-3946 [11]. The most promi- 
nent one has ~ 400 M® and is 1/10 of the SNR angular 
size. Thus, it should have a volume of (1 pc) 3 and therefore 
a column density of 20 /im of lead, that cannot produce a 
significant attentuation of the gamma rays. It has been ar- 
gued that the structure of this system of molecular clouds 
and the SNR shock waves are intimately connected [12]. 

Another feature of this SNR is that 7-rays up to 30 TeV 
or perhaps even more have been measured by HESS [13]. 
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Fig. 7. Various predictions for the muon flux per km 2 X yr expected 
from RX J 1713. 7-3946 and above 50 GeV. In blue, we quote the error 
deduced from the publications quoted in the text, in red, we quote 20% 
systematic error. In the horizontal axis, the year of the publication. 
The reasons of the changes are discussed in the text. From [8]. 

This is a pretty unique circumstance, in view of the limited 
area of the existing 7-ray telescopes that explored the TeV 
region but barely begun the exploration of the sky above 
few 10 TeV. The spectrum is non-trivial: it is not described 
by a simple power law, but rather by a broken-power-law 
or by a modified-exponential-cut as in Eq. 5 and with [14] 

7 = 1.79 ± 0.06 and E c = 3.7 ± 1 TeV (9) 

this does not seem to contradict the expectations, but the 
cutoff suggests that the cosmic rays with energies above 
150 TeV have already abandoned this SNR. 

5.1. Expectations within the hadronic hypothesis 

Thanks to HESS measurements, the upper bound on neu- 
trino signal is now quite precise. The most detailed recent 
calculation [2] deduced the and from these data, find- 
ing for a threshold of 1 TeV an induced muon flux of 

I^+f, = 2.4 ± 0.3 ± 0.5/km 2 yr (10) 

In other words, the HESS data, with the hadronic hypothe- 
sis, permit us a reliable evaluation of the expected fluxes-or 
more plausibly, the upper bounds-on high energy neutri- 
nos. It is not possible to lower much the detection threshold 
since the background due to atmospheric muons increases 
rapidly [4]: in fact, the flux of muon neutrinos decreases as 
~ E~ 3 ' 7 in the relevant region, as we see from Tab. 1. 

The expectations changed in the course of the time as 
shown in Fig. 7; let us recall the reasons. Ref. [15]: The 
first work on high energy neutrinos from RX J1713. 7-3946 
used a bold extrapolation of Cangaroo data, and adopted 
a simplified approach, subsequently improved. Ref. [16]: 
the inclusion of oscillations, Earth absorption, and livetime 
caused a factor of 4 decrease; moreover, the hypothesis of 
an early cutoff in the spectrum was remarked, estimating 
the maximum decrease in a factor of 3. Ref. [3]: the results 
of HESS, which first revealed a cutoff in the spectrum, per- 
mitted us to quantify the decrease in about a factor of 2. 
Ref. [2]: None of the subsequent theoretical and observa- 
tional improvements considered introduced large changes, 
and the prediction became stable. 

Two remarks are in order: (i) the last three predictions 
(or better, it is good to repeat it, upper bounds) agree 
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Fig. 8. The region around RX J 1713.7-3946 as seen by Fermi-LAT, 
superimposed with the location of the dense condensations named 
A,C,D seen by NANTEN-small circles marked by the red letters. 
The "Sources A,B,C '-larger circles marked by black letters-are the 
new point sources found by Fermi-LAT instead. From [17]. 

within errors; (ii) for a fair comparison, the threshold has 
been set to the common value of E t h = 50 GeV, but this 
value it is too low to permit background rejection. As shown 
in Tab. 1, this implies that the expected signal should be 
more than halved. 




Fig. 9. Mixed composition of the spectrum of RX 31713.7-39^6 as 
predicted by the Zirakashvili & Aharonian model. Figure from [19]. 

emission and a very efficient acceleration. From a theoreti- 
cal point of view, the spectrum could be leptonic (E~~< => 
E-y ( 7+1 )/ 2 ) or hadronic, with an energy-dependent pene- 
tration factor (E~~< =>■ E 7 7+1 / 2 ) 5 if primaries-electrons or 
protons, respectively-have 7 ~ 2 [12]. Moreover, the lack 
of thermal X-rays could imply that the emission is leptonic 
if the medium is uniform, or that non-uniformity plays a 
major role [18,12]. 

Since we expect progress from the GeV region in the close 
future, we believe that the right attitude is: Wait and see. 



5.3. Next steps? 



5.2. Fermi-LAT view at GeV and above 

However, the true point is whether the hadronic hypoth- 
esis holds true (i.e., whether the observed 7-rays come pre- 
dominantly from CR collisions) ; crucial tests of this hypoth- 
esis (that are of course relevant to the matter of predicting 
the high energy neutrino signal) include the observations 
of 7-rays in various regions of energy. 

Thus it is very interesting that, recently, Fermi-LAT re- 
leased their first results on RX J1713. 7-3946 [17], that re- 
vealed a rather complicated environment: 1) there is a wide 
source approximatively in the SNR location with spectrum 

I 7 oc E-^^ (11) 

that results from measurements above 5 GeV and from 
an upper bound on the spectrum at 0.5-5 GeV; 2) there 
are several point sources in the surrounding, including one 
sloping as E~ 2A5 , outshining the wide source at GeV 
energies; 3) finally, there is the diffuse background from the 
Milky Way. 

On the one hand, the location of the extended source does 
not disagree with the location of the dense condensations as 
shown in Fig. 8; however, what seems most crucial for the 
interpretation is the slope of the spectrum, that is different 
from the one given in Eq. 9 and also harder than theoreti- 
cally expected one, assuming that the observed emission is 
due to CR collisions. 

Let us pause here to recollect. From a purely observa- 
tional point of view, it seems very important to understand 
better the emission below 5 GeV, and one should note that 
a spectrum like E~ i:i is not yet excluded firmly-which 
could be reconciled more easily, perhaps, with hadronic 



The possibility that the spectrum of RX J1713. 7-3946 
is composite was discussed in a concrete theoretical model 
[19], see Fig. 9 (note however that the need to disentangle 
hadronic and leptonic emission was pointed out immedi- 
ately, see e.g., the statement of Aschenbach in [20]). This 
shifts the discussion on the need to quantify hadronic and 
leptonic 7 emission, rather than excluding one model in 
favor of the other one. Consider the possibility that the 
spectrum is composite. If we want to check thoroughly the 
SNR+MC paradigm, thus improving the expectations on 
high-energy neutrinos, we should not study only the spec- 
trum: we should also measure precisely the arrival direc- 
tions of the high-energy 7's, testing correlations with molec- 
ular clouds. This is especially important for 7's at 10 TeV 
and above, that are unlikely to be leptonic. Is it possible 
to use the existing data for this purpose? HESS has about 
500 events above 30 TeV, but 70% of them is expected to 
be cosmic ray background (i.e., noise) [13]. Then, in or- 
der to perform this check, we should wait for CTA or for 
an equivalent km 2 class 7-ray telescope, sensitive above 10 
TeV, and with good angular pointing. 

6. Other possible sources of high energy neutrinos 

Vela Junior: Above the bound necessary to have more 
than 1 /^/(km 2 xyr) [6], namely, 

7 7 (20 TeV) = (2 - 6) x 10~ 15 /TeV C m 2 s (12) 

there are 2 more young SNR, Vela Jr and Vela X, observed 
by HESS. The first one, also known as RX J0852-4622, 
is a young shell-type SNR, just as RX J1713. 7-3946, with 
angular size 2°. The specific parameters are however quite 
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unknown; the most common values of the estimated age 
are of 660-1400 yr, implying a distance of 0.26-0.50 kpc. 
Actually, this source is even more intense than RX J1713.7- 
3946 in 7-rays: 



/ 7 (20 TeV) = (1 - 3) x 10~ 14 /TeV cm 2 s 



(13) 



but 20 TeV is the last point presently measured by HESS. 
HESS spectrum oc E~ 2A , whereas Fermi-LAT's prelimi- 
nary measurement gives oc E~ 19 [21], which could indicate 
a spectral break. 

Cygnus Region: At 1.7 kpc from us, there is a large 
star forming region of 100, 000M Q mass in the direction 
of Cygnus constellation. This is remarkable target for cos- 
mic ray collisions, even if we do not know for sure which 
should be the associated source of cosmic rays. It includes 
several sources of TeV 7-rays, seen by Milagro [22], which 
could also radiate neutrinos. These arc potentially visible 
from IceCUBE location and two of them deserve a special 
mention: 

- MGRO 2019+37 is still unidentified and there is no visi- 
ble correlation with matter excess. ARGO & Veritas do 
not see it. Fitting 



I-y = 10~ n x £T 2 - 2 x e 



-y/E/E c 



/TeV cm 2 s 



(14) 



with E c — 45 TeV, it could imply up to 1.5 muon events 
per km 2 year above 1 TeV. 

MGRO 1908+06 is seen also by ARGO«Milagro>HESS; 
a pulsar has been found by Fermi-GLAST there. Using 



L = 2x 10~ n x E 



-2.3 



x e 



-y/E/E c 



/TeV cm 2 s (15) 



with E c — 30 TeV, it could imply up to 2.5 muon events 

per km 2 year above 1 TeV. 
A few remarks: 1) The cutoffs E c have been introduced to 
account for the CASA-MIA bounds at 100 TeV. 2) There 
is another source, MGRO 2032+41, but weaker. 3) Let us 
repeat that these objects are weaker theoretical cases than 
SNR's, but at least, they include a large target material. 
Fermi Bubbles: Possible cases of diffuse/wide sources 
could be given by the recently discovered 'Fermi bubbles'. 
Are they a reservoire of galactic cosmic rays? If so, they 
could be also promising neutrino sources [23] . By extrapo- 
lating the flux to 



I 7 {E) = Q IQ^e-^^/E 2 / TeV cm 2 s 



with E c =100 TeV (meaning a cut at 1 PeV in CR spec- 
trum) and using an angular size of SI — 0.2 sr m n x (15°) 2 , 
we get a signal of about 100 muons a year for 1 km 2 detec- 
tor area. This could be observable in Km3NET as a wide 
neutrino source over the background. 

7. Summary and outlook 

The promising galactic neutrino sources are intimately 
tied to 7's above 10 TeV. The main points regarding the 



predictions of the high-energy neutrino fluxes from galactic 
sources are simply two: 

- There are only few bright 7-ray sources, that correspond 
to detectable neutrino sources. 

- For these sources, the expected neutrino signal is mod- 
est, even assuming that the 7-ray emission is fully con- 
tributed by hadronic processes. 

The studies of 7-rays from RX J1713. 7-3946 prove that it is 
possible to proceed toward reliable expectations for high- 
energy vs. The future years will offer us good occasions of 
progress in this sense: 

(i) Multiwavelength observations and theory will help us 
to understand better the young SNR's and to quantify 
their hadronic and leptonic emissions. 

(ii) Sub-degree pointing with very high energy 7-rays 
could permit us to further test the paradigm of 
SNR+MC association for RX J1713.7-3946. 

These considerations emphasize the relevance of certain 
lines of progress, and suggest to increase synergies and col- 
laboration between neutrino and gamma ray search. 

Appendix A. Neutrino oscillations 

Averaged oscillation probabilities: Three flavor oscilla- 
tions are well-understood and relevant. For the transparent 
sources we are discussing, the simplest regime [24] applies 
and the relevant probabilities are just constant: 



(A.l) 



In order to provide explicit expressions of the averaged 
probabilities, we adopt the standard parameterization of 
the leptonic mixing matrix to solve the unitarity contraints. 
Then, defining 



e = sin t/13 

a = sin 813 cos 5 sin 2# 12 sin 26* 23 
/3 = (1 -e) cos2<9 12 /2 + 7 
7 = (1 + e) cos 26»i2 cos 26> 23 /2 



(A.2) 
(A.3) 
(A.4) 
(A.5) 



we obtain relatively simple expressions for the averaged 
probabilities: 



(16) Pee = (1 - e) 2 1 - - sin 2 20 12 + e 



1 



efi - 



P — 



[cos 2 6*23 sin 2 26»i2 + 4e sin 2 6*23 x 



sin 2 2012 , 
1 : I + a. cos 26 ) 12 



1 3 

3 + 2 



(1 - e) sin 2 023 



+ -(a-py 



(A.6) 

(A.7) 
(A.8) 



P — 

1 fin 



[(1 + e) 2 + (1 - e) 2 (3 sin 2 26 23 - sin 2 20 12 )] + 



(A.9) 



Fig. A.l. Dependence of the effective oscillation probability on the 
CP-violating phase and other oscillation parameters; see the text for 
the meaning of the various lines. 

One can check explicitly that the diagonal probabilities 
obey 1 > Pu > 1/3 and the out-of-diagonal ones 1/2 > 
Pli 1 > 0. The expressions of the remaining probabilities 
P eT and P TT are obtained from these expressions, mapping 
5 — >• S + 7r and 6*23 —> tt/2 — #23. Note that these formulae 
depend on e, on the 'doubled angles' 2#i 2 and 28 2 3, and on 
the combination sin #13 cos S. 

The muon neutrino/antineutrino flux after oscillation 
Fp, as a function of the fluxes before oscillation, F® , is: 



P — p pO , p pO — p jpO 



(A.10) 



Two important remarks on P eff are in order: 

1) P eff depends on the energy through the flux ratio F®/F®. 
A typical value for CR collisions is 1/2, for when we consider 
an almost equal amount of charged pions tt 1 ^ , 1 electron 
(anti)neutrino each 2 muon (anti)neutrinos are produced, 
all with similar energies (in [2] one finds more precise state- 
ments for power law distributions). However, other cases 
are possible at least at a speculative level: e.g., if muons 
are absorbed before decaying (which could happen, e.g., 
in microquasars) the initial beam is purely made of muon 
neutrinos, F®/F® — 0; or, if the neutrino originate from 
neutron decays, only electron antineutrinos are produced 

2) The probability P cS depends also, to some extent, on the 
most poorly known oscillation-parameters. Fig. A.l shows 
the variation of P e s with the CP-violating phase 6, where 
#12 = 34° is well-known and thus fixed, whereas the other 
angles are moved in the wide and conservative ranges #23 = 
(40° -> 50°) and 13 = (6° -> 12°)-wiggly yellow band- 
and we set F® /F® = 1/2. For comparison, we show also the 
curves obtained by setting 623 = 50° and #13 = 6° (upper, 
red) and 6*23 = 40° and 6» 13 = 12° (lower, blue); the lower 
part, grey region, is forbidden. 613 is a relevant variable, 
thus future precise measurements will impact somewhat on 
oscillation probabilities. 



Numerical values: The values of Pw change a bit with the 
new best fit values of 6*23 = 40° and #13 = 9° of [25] (assum- 
ing the new reactor fluxes) along with [26] 5 — —90°, that 
replace the previous reference (or extremal) values #23 = 
45° and #13 = 0°, now disfavored; 6*12 changed slightly 
35° — > 34° . The new values imply that the symmetric ma- 
trix with elements Pw changes as follows; 



P 

0.40 



Fig. A. 2. Values of P^^ and P ell , compatible at la, 90% andla K, 95% 
with the present information on the mixing angles, as described in 
Eg. A. 12. 
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The changes are not very large, though. 

The present ranges of parameters, approximatively 95%, 
are 6 12 = 32° - 36°, 6 13 = 4° - 12° and (9 23 = 37° - 51°, 
compare [25] and [26]. The precise values of the best fit 
and of the ranges depend on the details of the analysis. 3 
For illustration, we adopt in the following the comparably 
larger value of 613 and smaller central value of #23 of [25], 
and take into account also the existing (feeble) indication 
for S < [25,26], assuming: 



6» 12 = 34°±1°, sin 2 13 = 0.025 ± 0.007, 
9 23 = 42° ±3.5°, 



2 4 



(A.12) 



that allow us to perform simple statistical analyses. More 
precisely, we can consider these parameters as independent 
Gaussian variables, and constrain the averaged probabili- 
ties as described in [20] . 

Now, let us consider the allowed ranges for the two 
most important averaged probabilities, P eu and P U n- From 
Eq. A.12, we derive at la the values P PAl = 0.37l 03 and 



P. 



0.25 



+0.04 
-0.05- 



Fig. A. 2 shows that these two quantities 



are strongly correlated; thus, it is not immediate to derive 
the expected range of P e g. This can be done however prop- 
agating the error directly; the result is shown in Tab. A.l. 
For the important case F® /F° = 1/2, we find the 2a range 
P efi = 0.32 — 0.36, which is wider but does not disagree 
severely with the result in [2], P cff = 0.33 — 0.35. 

All in all, the new results may eventually lead to a small 
increase of P e g, and thus of the signal: about 5%. In view of 
the larger particle physics uncertainties, and much larger 
astrophysical ones, we maintain our position, expressed in 
2004 [16], that the small uncertainties in the oscillation pa- 
rameters imply that "there is little hope to learn anything 



3 In particular, there is an important dependence on the (theoretical 
and experimental) value of the flux from the reactors, discussed 
in [26], that is another reason why the measurements of Double 
CHOOZ, Reno and DAYA Bay will be important to proceed further. 



useful on 3 flavor oscillations" studying high energy neu- 
trinos from cosmic sources. Different opinions have been 
however expressed several times in the literature, and au- 
thoritatively, in [27] ; our criticisms to this specific proposal 
can be found in [20], footnote 1. 

Expansion in small parameters: It is also possible to ex- 
pand the averaged probabilities P«< in sin 9\z and cos(2# 2 3) 
that are small parameters and, due to the present uncer- 
tainties, vary in a similar range; in fact, sin(10°) = cos(2 x 
40°) = -cos(2 x 50°) ~ 1/6. Thus, we expand the aver- 
aged probabilities around the point 6*i 3 = and 9 2 z = 45° 
(namely, where the two small parameters are zero) neglect- 
ing cubic terms and finding: 



Pee = 1 - 2 - 2W > *V 



X „ 1 X 

Pe l l=^+V + W , P^ = -- -- y-W + Z 



(A.13) 



X 



P eT = -- y + w, P T 



x 



y — w + z 



We introduced for convenience various quantities of differ- 
ent size, 



order 1: X '■■ 
linear: y ■ 



:sin 2 (26> 12 ) 
1 



(A.14) 



quadratic: Z \ 

(l-s/4) 



sin #13 cosb\/ x(\ — x)+cos(26>23)x (A. 15) 

(A.16) 

sin 2 6»i3 (l+cos(2£)a;/2) + cos 2 (26> 23 ) x 

(A.17) 



quadratic: W = sin #13 (1 — x/2) 
1 
2 

sin#i 3 cos(2#23) cos6y/x(l — x) 



the linear expansion has been obtained already in [16] . Sev- 
eral remarks are in order: 

(i) By summing x, y, w or z in any row or column of the 
symmetric matrix with elements P u > one gets zero; 
this is a consequence of unitarity. 

(ii) The unique linear term y involves a weighted sum 
of sin6>i3 and cos (2023); this implies that there is a 
strong correlation between the linear variation of P e ^ 
and P w , while P ee and P MT do not vary at this or- 
der. Note that due to the coefficients, the impact of 
cos(26*23) is larger than the one of sin 6*13. 

(iii) The relative size of y and z depends on 5; the latter 
can be larger than the former (when <5 ~ 7r) even if 
these two terms are formally of different orders. 

These expressions can be useful for some purposes, e.g., to 
understand the correlations among the probabilities, how- 
ever they are not much simpler than the full expressions. 
Moreover, the a priori limits 1 > P« > 1/3 and 1/2 > 
P u , > (when I ^ £') and the specific values of the param- 
eters of oscillations imply that P eff is close to the physical 
boundary [2]; in other terms, the non-linear terms have a 
certain role, as consistent with the previous last considera- 
tions. For these reasons, we decided to use the full expres- 
sions of Pw for the above analysis. 
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Table A.l 

Value of the effective probability of oscillation, Eq. A. 10, and of its 
error, for various values of the flavor ratio F® /F®. 
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